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Abstract 

Performance of two dimensional foil sections with traditional, blunt and cupped trailing edge 
geometries are experimentally and computationally evaluated. Traditional foils show less than 
1% error in lift slope and less than 0.1 degree error in predicted angle of attack. Foils which 
include trailing edge separation due to a trailing edge cup show up to 30% error in lift slope 
and 1.0 degree error in predicted angle of attack. Foils which include trailing edge separation 
due to bluntness show good correlation in lift slop (2% error) but still show up to 0.9 degree 
error in predicted angle of attack. 

Experimental and numerical evaluations are conducted in order to assess whether the differ- 
ences are caused by the experimental or computational fluid dynamics. All known experimental 
uncertainties are exhausted without explaining the differences between predicted and measured 
lift. 
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Title: Professor of Naval Architecture 
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Chapter 1 



Introduction 



1.1 Motivation 

Military and civilian ships are demanding higher speeds and improved cavitation performance. 
These requirements are pushing current propulsion system designers to rapidly improve existing 
systems. One of the directions these advancements are pushing designers is toward novel new 
section designs in propeller blades. 

Based upon analyses accomplished to date, it is not clear computational fluid dynamics 
codes are accurate for the advanced section designs considered. The following is the abstract 
from a thesis conducted by Jorgen Jorde [l] in December 1995: 

“This thesis covers a comparison between experimental and numerical results of 
a foil with a slightly cupped trailing edge, and a foil with a more distinctly cupped 
trailing edge. 

Experiments have been conducted to find lift and drag coefficients, as well as 
boundary layer parameters at the trailing edges of both foils. The experiments 
indicate that the numerical tools available are not capable of accurately predicting 
the lift of these foils. 

Analysis of the results further indicate that experiments where there is a great 
wall/foil interaction are difficult to interpret until numerical tools are refined to the 
extent that they may be the link to the performance in unbounded flow.” 
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Figure 1-1: Range of Foil Geometries Considered 

The numerical tools have now been refined to the extent that they can analyze the behavior 
of multiple foil elements and simulate the effects of the tunnel walls. The present thesis 
attempts to answer the questions raised by Jorde by applying the advancements in numerical 
tools over the last five years and improving upon the experimental techniques and analysis. 



1.2 Goal 

The goal of this thesis is to assess the capabilities of experimentation and computational fluid 
dynamics in predicting performance of challenging foil section designs. 



1.3 Method 

This thesis initially focuses on a single foil, at a single angle of attack, at a single Reynolds 
Number and seeks to assess whether there is agreement between experimental and computa- 
tional analyses of this foil under this specific set of conditions. The foil with the highly cupped 
trailing edge shown in Figure 1-1 is the initial area of focus and is the same foil used by Jorde 
[l] analyzed in 1995. Questions that arise from the analysis of this foil drives the thesis toward 
a wider consideration of the foil geometries depicted in Figure 1-1. 
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1.3.1 Cupped Foil Experimental Analysis 

This section improves upon existing experimental methods and applies them to establish a 
confidence interval for the performance of the foil being assessed. The initial assessment of 
the experimental methods available call into question the ability to accurately measure fluid 
dynamic drag. Methods are developed for determining pressure in the wake which allows re- 
finement of traditional methods as well as a new thermodynamic analysis to determine drag. 
The methods for calculating drag are applied to a control case of a two dimensional cylinder to 
assess their performance. These methods are then applied to this foil to establish a confidence 
interval for fluid dynamic drag. The initial assessment of the experimental methods for cal- 
culating lift seem to indicate that existing methods are acceptable. An extensive data field is 
taken to allow a convergence study, error analysis and confidence interval for the fluid dynamic 
lift. 

1.3.2 Cupped Foil Numerical Analysis 

This section applies the most recent updates of numerical procedures to establish a confidence 
interval for the performance of the same foil in computational fluid dynamics. This section 
shows good agreement between computational fluid dynamics methods, but the predicted lift 
slope is 30% larger than the measured lift slope and there is a one degree offset in angle of 
attack. This discrepancy between experiment and computation is so large that it brings into 
question the experimental results as well as the computational results. 

1.3.3 Analysis of Numerical and Experimental Solution Incompatibility 

This section of the thesis seeks to explain why the experiment and predictions are so far out of 
agreement on this foil. Several avenues that did not bear fruit in explaining the differences are 
discussed. At the conclusion of the analyses on the very difficult highly cupped foil both the 
experiment and computations appear valid, but still show very different answers. 
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1.3.4 Analysis of Other Foil Sections 



This section of the thesis then seeks to see whether the performance of other foils in our water 
tunnel can be predicted. The live foil geometries shown in Figure 1-1 are evaluated. In this 
section it is shown that there is very good agreement between the computational and experi- 
mental lluid dynamics for the case where the foil takes the traditional approach of unloading the 
foil over the last 10% of the chord length and ending at a sharp point. It is shown that when 
the foil ends at a blunt trailing edge or when it is not unloaded at the trailing edge agreement 
is not achieved. 

1.3.5 Conclusions 

In the conclusions the assertion is made that CFD performs well for foils that are unloaded 
near the trailing edge, but shows significant errors when separation occurs under asymmetric 
loading. Some speculative reasons for these inaccuracies are proposed. 
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Chapter 2 



Cupped Foil Experimental Analysis 



The experimental analysis conducted by Jorgen Jorde [1] is extensive. He conducted a sig- 
nificant evaluation of the gross flow field velocities and boundary layer measurements. The 
experimentation and analysis conducted here seeks to extend and improve upon his efforts to 
allow more accurate determination of drag coefficient as well as error analysis for both lift and 
drag. 

2.1 Experimental Overview 

The MIT Variable Pressure Water Tunnel (Figure 2-1) is a closed loop system. A variable 
speed impeller forces flow through the loop at speeds up to 9 m/sec. 

The test section is twenty inches square and features removable clear windows on all four 
sides. In the test section, a three watt argon-ion laser doppler velocimeter is used to measure 
local vertical and horizontal velocities. The laser measurement volume is less than 0.07 mm 
across and 0.7 mm deep. Mean velocities are repeatable to within 0.002 m/sec with the sample 
sizes taken in this experiment. Velocities are accurate to within 0.6%. A Kiel Probe is 
used to measure differential stagnation pressure from the front of the foil to the wake. The 
measurement diameter for the Kiel Probe is approximately 2 mm. 

Figure 2-2 is a depiction of all data taken in the initial phase of this experiment. The 
chord length of the foil is 455 mm. Fluid velocities are taken at a path length interval of 3 mm 
along all of the lines except for the vertical lines in the wake of the foil. The vertical lines in 
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Figure 2-1: MIT Variable Pressure Water Tunnel 

the wake of the foil are actually 5 vertical lines of data taken at a 1 mm spacing which allow 
for calculating flow gradients in the longitudinal direction. The foil angle of attack is fixed 
at 0.35±0.05 degrees. The horizontal tunnel velocity is maintained constant throughout the 
experiment at 7.003 m/sec. The velocity is held constant to within 0.3% and is accurate to 
within 0.6%. 

2.2 Pressure Measurement 

The pressure throughout the majority of the fluid domain is determined based on Bernouli’s 
Equation [2]. In the wake there has been significant fluid flow energy converted to internal 
energy such that Bernouli's Equation is not valid. In this region two methods of pressure 
measurement are used in this experiment. The first method, already discussed is the Kiel 
Probe. The second method is to apply measured velocity gradients to the Navier Stokes 
Equations and calculate the pressure variation. 
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Figure 2-2: Test Section Layout and Data Summary 



2.2.1 Deri vat io n 

Beginning with the vertical component of the Navier Stokes Equation, 



( 2 - 

where 

u = fluid velocity in x direction 

v = fluid velocity in y direction 

w = fluid velocity in z direction 

p = static pressure 

p = fluid density 

ft = fluid viscosity 



Navier Stokes Equation 

civ dv dv 1 dp 

U— + v— + w— = — + pSJ V 

dx ay dz p dy 



The flow domain considered for this application is approximated as two dimensional so the 
term is negligible. The low viscosity of the fluid integrated over this small length causes 
the viscous (p\7 2 v) term to also be negligible. Each of the velocities is considered with a mean 
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and fluctuating component which yields Equation 2.1 from Equation 2.2. 



1 dp 

pdy 



_ dv _ dv , dv , dv 

u - — f U— h U-~ — I- v— 

dx dy dx ay 



where 



(2.2) 



u — mean x direction velocity 
v = mean y direction velocity 
u = fluctuating x direction velocity 
v = fluctuating y direction velocity 

The fluid velocity gradients in the x direction are determined using a least squares fit from 
five points at an interval to give best gradient resolution. In order to transform this equation 
into a form reliant on horizontal gradients and into a form reliant on physically measured 
quantities, a few substitutions are made. Equation 2.3 is simply the continuity equation in two 
dimensions. 



du dv 
dx dy 



Equation 2.4 is simply an identity for the derivative of a product. 



(2.3) 



, dv 



dx 





(2.4) 



Applying these equations and taking the time average gives Equation 2.5 which is the form 



that is directly applied to the measured data. 



1 dp 
pdy 



_ dv _ du d — d ts- 
= u - — h v — — b — uv + — v z 
dx dx dx dy 



(2.5) 
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Figure 2-3: Pressure Variation by Component 

The laser doppler velocimeter used in this experiment simultaneously samples the velocity 
of a particle in the horizontal and vertical directions. This feature allows a statistical analysis 
of uv and v 2 at each point in space. 

2.2.2 Results 

Appendix B shows the raw data taken in this experiment. The ^uv term is found to be two 
orders of magnitude smaller than the other components. The effect of the other individual 
terms shown in Equation 2.5 on pressure is plotted in Figure 2-3. The has very little 
impact on the overall problem which could have been anticipated by the fact that both v and 
2 “ are small quantities. The term provides the appropriate slope for the two sides of the 
boundary layer to return to the expected Bernouli Pressure. The — v 2 term gives the pressure 
dent in the center of the wake that is caused by the turbulent viscous effects. 

Figure 2-4 shows the results of plotting stagnation pressure measured with a Kiel Probe 
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Calculated from Speed Data 

o Measured with Ki el Probe 

Figure 2-4: Stagnation Pressure Coefficient vs Vertical Wake Position 

and stagnation pressure inferred from the Navier Stokes Equations as shown above. 

As a brief comment to experimentalists, obtaining accurate flow gradients in the wake region 
is not a trivial matter. In order to obtain sufficiently repeatable vertical velocities to accurately 
distinguish gradient information it was necessary to take 3000 velocity samples at each point in 
space. On top of this, it was necessary to take five separate sample points at an appropriate 
spacing in the horizontal direction to obtain the slope using a least squares fit. In short, 
this means 15,000 velocity samples to obtain flow gradients at a single point in space. The 
gradients through the wake need to be resolved at a maximum spacing of 1 nun. With a wake 
thickness of about 40 mm this translates to 600,000 velocity samples to obtain the pressure 
information across the wake at one location. Figure 2-5 shows the pressure determined from 
velocity magnitudes and gradients at three longitudinal positions in the wake. 
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Figure 2-5: Stagnation and Static Pressure in the Foil Wake at Three Longitudinal Positions 

2.3 Drag Measurement 
2.3.1 Overview 

Experimental measurement of profile drag is given considerable attention in wind tunnel testing 
and water tunnel testing. Measurement of profile drag requires knowledge of pressure and 
velocity. Until recently, velocity was inferred from pressure measurements because direct 
measurement of fluid velocity was not available. Sclilichting [2] summarizes various methods 
of calculating drag on a foil based strictly on pressure measurements in the wake of the foil. 
These methods have been demonstrated to be accurate when taking measurements in the wake 
of the foil at least 0.05 chord length behind the foil. A potential weakness of these methods 
explored in this thesis is inaccuracy dealing with unsteadiness in a wake. 

In recent years laser doppler anemometry has given experimentalists the ability to directly 
measure velocity in a fluid domain. This capability provides great insight into the nature 
of the flow in the wake and allows accurate measurements of fluid velocity and turbulence 
intensity. Kinnas [3] developed a method for inferring fluid dynamic drag based solely on 
velocity information. This method assumes the pressure across the wake can be determined 
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Figure 2-6: Velocity Definitions for Analyses 

by linear interpolation from one edge to the other. The difference in static pressure between 
the method of Kinnas and the measured values shown in Figure 2-5 equates to an 8% error in 
drag measurement. The bigger problem with this method for this particular application is the 
effect of sidewall boundary layer growth on the overall measurement. Appendix A summarizes 
the effect of sidewall boundary layer growth on a control volume analysis for drag. It has been 
the conclusion of Kimball [4], Jorde [l] and this study that the bounding box method is not 
well suited to water tunnel testing. 

This study adds to these methods another based on thermodynamics. The theory is covered 
and then a very challenging unsteady problem of measuring drag on a cylinder is tested. The 
thermodynamic method provides results within 2.5% of the direct measurement of the cylinder 
drag while all of the other methods are 13% or more in error. Three methods are then applied 
to the foil of interest and the average value of each method is found to agree within 1% of the 
other methods. The reason the methods agree well in the case of the foil, but not in the case 
of the cylinder is that only the thermodynamic method takes into account unsteadiness in the 
wake. The cylinder has significant unsteadiness in the wake and the foil has relatively little. 



20 



2.3.2 Experimental Drag Measurement Theory 



Momentum Methods 



Equation 2.6 can be used to determine fluid dynamic drag on the foil if the fluid velocities 
and pressures are known all the way around a closed contour that surrounds the body. One 
example of such a contour is the flow domain that stretches from the top to the bottom of the 
tunnel in front and in back of the foil. This control volume is the one used in the method of 
Betz and the method of Jones, which are summarized by Schlichting [2]. Another possibility 
would be to simply measure the fluid velocity all the way around the foil, and this control 
volume is the one used in the method of Kinnas [3]. 

Momentum Contour Integration (Kinnas) 

D = -f>§ [uvdx + (u 2 + u 2 ms ) dy\ + f pdy 
u=horizontal velocity 

( 2 - 6 ) 

v ~ vertical velocity 
p— static pressure 

Analyzing the control volume that stretches from the top to the bottom of the tunnel is 
complicated by the fact that velocities are not known all the way to the edge of the tunnel. In 
order to analyze this control volume it is necessary to transform Equation 2.6 into a differential 
form that tends to zero outside of the wake. In Betz formulation, a false source is used to raise 
the fluid velocity u2 to u2’ to provide a useful flow field with a constant stagnation pressure. 
He then uses the drag associated with the false source as a correction to the change in total 
pressure portion of the formulation. Equation 2.7 is used to calculate fluid dynamic drag based 
on the method of Betz. 

Betz’ Formulation 

D = J'Zc \ai - .92] dy + § [(«2 - u' 2 ) K + u 2 - 2 C/00)] dy 

g = p + ef- (2.7) 

£/,x>=free stream velocity 

Jones formulation assumes the fluid will travel from 2 to 3 with no further losses. Point 3 
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is sufficiently far downstream to allow the static pressure to equalize across the wake. He 
uses continuity along a streamtube to project conditions at 2 back to point 3 where the static 
pressure has returned to free stream conditions and is no longer contributing to drag. He 
then simply evaluates the momentum deficit at point 3. Equation 2.8 is used to calculate fluid 
dynamic drag based on the method of Jones. 



Jones’ Formulation 

/ oo 

[u 2 (U , o - u 3 )} dy (2.8) 

OO 

Thermodynamic Method 

The first law of thermodynamics states that the rate at which work is done on a control volume 
is equal to the net energy flux out of the control volume. Consider the control volume between 
1 and 2 in Figure 2-6. The force on the foil is stationary relative to a fixed coordinate system; 
therefore, no net work is done on the control volume. Based on the first law of thermodynamics 
there can be no net change in fluid energy from 1 to 2. There is a reduction in flow energy 
and kinetic energy that is being converted to fluid internal energy. This process is irreversible. 
Equation 2.9 shows the change in fluid internal energy from 1 to 2. 

^2 — (^^oo T Poo^j — + £oo — 9oo — 92 + ^oo (2.9) 

E ~ Fluid Internal Energy (2-10) 

In order to evaluate the system using the second law of thermodynamics it is desired to project 
the change in internal energy to 3. Jones' method was to assume no further losses, but with 
flow separation unsteadiness in the wake must also be considered. It is assumed that the energy 
associated with the unsteadiness in the wake will be converted to internal energy via viscous 
dissipation. Calculation of energy associated with the unsteadiness is taken from Schlichting 
[2]. The fluid internal energy at 3 is then given by Equation 2.11. 



“ *?oo 92 



i (“ 2 



rms 2 



+ V, 



rms 2 



) + E, 



J oo 



( 2 . 11 ) 
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Equation 2.12 is the basic relation for calculating a change in fluid entropy based on a 
change in fluid internal energy. It is a direct result of the second law of thermodynamics [5]. 



TdS 



P 

v 

S 



dE + Pdv 
pressure 
change in volume 
entropy 



( 2 . 12 ) 



Water is incrompressible so there is no change in volume so Equation 2.12 reduces to 2.13. 



^3 *§oo — 



£3 — Eo 

n~ 



(2.13) 



The ‘lost work rate” associated with this entropy generation is given by Equation 2.14. 



W lost— j [To (S 3 - Soo) u] cly 

The drag on the hydrofoil based on this analysis is given by Equation 2.15. 

tj Wiost 

~ U 



(2.14) 



(2.15) 



2.3.3 Results 

In order to assess the quality of the various methods in measuring fluid dynamic drag in an 
unsteady wake a cylinder is evaluated as a test case. The actual drag on the cylinder is 
determined by measuring the pressure on the surface of the cylinder through a very small hole 
and rotating the cylinder. This method of drag measurement assumes small tangential viscous 
forces on the cylinder surface. The static pressure is integrated around the cylinder to obtain 
the drag force. The experimental method and results are recorded by Prieto [6]. The direct 
measurement and each of the other method solutions are summarized in Table 2.1. Each of 



the methods show significant error with the exception of the Thermodynamic method which is 
within 2.5% of the direct measurement. 



Technique 


Drag Coefficient 


Direct Measurement with Pressure Integration 


0.77 


Betz 


0.67 


Jones 


0.67 


Kinnas 


0.49 


Thermodynamic 


0.75 



Table 2.1: Summary of Method Results for Cylinder Drag Coefficient 



The fluid dynamic drag is also measured for the foil. The drag is calculated at 3 different 
longitudinal positions; the results are shown in Figure 2-7. The method of Kinnas, Betz and 
Jones were all applied using the calculated pressure discussed above and the measured velocities. 
As discussed in Appendix A the method of Kinnas diverges as the control volume size increases 
due to three dimensional effects. This accounts for the significant deviation between this 
method of measuring drag and the other three. The methods of Betz and Jones both use 
the assumption of convection without further losses; therefore, it makes sense that their results 
would agree. The thermodynamic method provides results very similar to the methods of Betz 
and Jones. A simple average for each of the three valid methods shows remarkably consistent 
results as shown in Table 2.2. The conclusion is that the drag coefficient for this foil is 0.0097 
with an error margin of 2.5% to encompass all experimental results and stay consistent with 
the errors measured in the cylinder experiment. 



Technique 


Drag Coefficient 


Betz 


0.00969 


Jones 


0.00973 


Ther mod}uiamic 


0.00973 



Table 2.2: Summary of Method Results for Foil Drag Coefficient 



2.4 Lift 

Equation 2.16 is derived in many textbooks such as Newman [7]. 
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Figure 2-7: Summary of Drag Calculation Methods for the Cupped Foil 



Momentum Contour Integration 

L = ~P § [uvdy + ( v 2 + vims) d x \ + § pdx 
u=horizontal velocity 

(2.16) 

v=vertical velocity 
p=static pressure 

Figure 2-8 shows the foil in the tunnel with the various control volumes analyzed for lift 
coefficient by applying Equation 2.16. Appendix B summarizes the data taken along these 
contours. Based on these results the lift coefficient is determined to be 0.462db0.01. 



25 




C l =0.002 


0.004 


C l =-0.005 


0.005 


0.000 



Figure 2-8: Summary of Results for Measured Lift Coefficient 
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Chapter 3 



Numerical Predictions of Cupped 
Foil Performance 



Two computational fluid dynamics codes are used to evaluate the performance of the cupped 
trailing edge hydrofoil for which experimental evaluation has been conducted at the MIT vari- 
able pressure water tunnel. Both codes used model both the foil and the tunnel wall in two 
dimensions. The two flow codes utilized are IFLOW [8] and MSES [9]. 

3.1 IFLOW 

IFLOW is a Reynolds Averaged Navier Stokes code developed by Sung [8]. IFLOW is capable 
of solving the two and three-dimensional incompressible steady and unsteady RANS equations 
and uses a non-linear k - u turbulence model for closure. The turbulence model is the standard 
k - uj model developed by Wilcox [10] coupled with a nonlinear Reynolds stress model. 

The cupped foil is analyzed in an infinite fluid domain and also with the tunnel walls. The 
foil performance is verified to be converged by iteration steps and mesh density. Table 3.1 
shows the results of this analysis. 



Boundary Condition 


Lift Coefficient 


Drag Coefficient 


Infinite Fluid 


0.525 


0.011 


With Viscous Tunnel Walls 


0.625 


0.012 



Table 3.1: Iflow Numerical Results for Cupped Foil 
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3.2 MSES 



MSES uses a boundary layer model to predict boundary layer growth and separation along the 
foil. The inviscid Euler formulation and integral viscous formulation are coupled through the 
displacement thickness and solved simultaneously with a global Newton method. The code is 
developed and adapted for use with blunt trailing edges by Mark Drela [9]. Three methods of 
analysis are used in MSES. 

1. Foil performance is measured in an infinite fluid domain. 

2. A built in feature for wind tunnel testing is used to provide the effect of inviscid tunnel 
walls on the flow regime. 

3. The ability to analyze multiple airfoil elements is utilized to put a long thin plate above 
and below the foil being analyzed. This long thin plate develops a boundary layer 
similar to the tunnel walls and therefore provides a more accurate depiction of the actual 
experimental setup. The results determined using this method are 5% lower than the 
results without top and bottom boundary layers. All results in the following sections are 
reported based on analyses using this method. 



Boundary Condition 


Lift Coefficient 


Drag Coefficient 


Infinite Fluid 


0.545 


0.009 


With Inviscid Tunnel Walls 


0.623 


0.010 


With Viscous Tunnel Walls 


0.590 


0.010 



Table 3.2: Mses Numerical Results for Cupped Foil 



3.3 Conclusions 

The computational analyses conducted here are about 6% different from one another and about 
and about 30% off of the experimental measurements. This represents an angle of attack offset 
(assuming experimental lift slope) of about 0.2 degrees between the two codes and about 1 
degree between the experiment and the codes. The results are outside the margin of error 
in the experiment. The question at hand is whether the margins of error in the experiment 
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have not been appropriately considered or whether the computational fluid dynamics codes are 
inaccurate. 



29 



Chapter 4 



Analysis of Numerical and 
Experimental Solution 
Incompatibility 



The cupped foil geometry shows consistently poor correlation between experiment and compu- 
tations over a wide range of angle of attack and Reynolds Number. Figure 4-1 is a summary of 
the lift measurements taken on this foil. As shown in the diagram, both the computational and 
numerical calculations have been performed by a range of a people, over a significant period of 
time, using different equipment and different algorithms, but the end result is a consistent lift 
prediction that exceeds measured lift prediction by 0.1 to 0.2 depending on angle of attack. 

This chapter poses a series of questions that have been addressed in an attempt to reconcile 
the differences between computational fluid dynamics results and experimental results. Each 
section includes the parameter or concept in question as well as how this question has been 
addressed. Several of the sections refer to particle image velocimetry data that was taken to 
augment the laser doppler velocimetry data. Appendix C covers the experimental procedures 
and data collected via particle image velocimetry for this experiment. 
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- 0.5 - 0.3 - 0.1 0.1 0.3 0.5 0.7 0.9 

Angle of Attack 



Figure 4-1: Lift Coefficient vs Angle of Attack for Cupped Foil 

4.1 Angle of Attack 

The first question that naturally arises is whether or not the same angle of attack is being used 
in the flow codes and in the experiment. Four independent methods of measuring the angle of 
attack are used. 

1. The vertical position of the leading arid trailing edges are located using the horizontal 
beam of the laser doppler velocimeter. 

2. The position of various points on the surface of the foil are measured with the laser and 
these points are fitted to the foil offsets rotated at a specific angle of attack. This method 
allowed an error analysis of the foil angle of attack by rotating beyond the fit to the point 
where it is clear that there is no longer agreement. Figure 4-2 shows the foil at 0.35 
degrees along with the offsets measured. It was found that if the foil was rotated to 0.3 
or 0.4 degrees it is clear that the offsets no longer agree. Based on the second method 
the foil angle of attack is determined to be 0.35±0.05. 

3. A mechanical pointer is connected to a traverse to measure the position of the leading 
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design foil offsets at 0.35 degree angle of attack 

OOO foil offsets measured with the laser 



Figure 4-2: Laser Varification of Foil Angle of Attack 

and trailing edges of the foil relative to the tunnel horizontal axis. 

4. Particle image velocimetry pictures are combined to show the position of the surface of 
the foil along its length. 

In summary, the angle of attack is confident within 0.05 degrees. From Figure 4-1 it is clear 
that the error in angle of attack is insufficient to be a significant contributor to the discrepancy 
between experimental and computational results. 



4.2 Foil Geometry 

The next question that naturally arises is whether or not the foil being modeled in the compu- 
tational fluid dynamics codes is the same foil analyzed in the experiment. The foil geometry 
is verified as follows: 

1. The foil offsets from the computational fluid dynamics codes are scaled to the experimental 
foil size. 

2. A computer controlled waterjet cutter is used to cut out the foil shape in aluminum. 

3. The foil shape is slipped over the foil and inspected for gaps at the edges 

The conclusion of this process is that the actual foil shape and evaluated foil shapes are the 
same. 
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Figure 4-3: Predicted Verses Measured Edge Velocity 

4.3 Data Processing Errors 

Another issue to consider is whether or not the lift coefficient for the experiment is calculated 
accurately. This concern is alleviated in three ways. 

1. The flow velocity information from I FLOW is used to set up a mock experiment. The 
fluid velocity at all of the points measured in the experiment are extracted and used to 
calculate lift using the momentum integration method used in the experiment. The result 
is that the mock experiment agrees with the IFLOW lift output to the third significant 
digit. This shows that in the case of perfect data our algorithm performs well within our 
predicted error margin. 

2. Figure 4-3 shows the edge velocities on the suction side of the foil verses chord position. 
Both the computational and experimental velocities are normalized by velocity measured 
upstream at a geometrically identical location. It is clear that both computational fluid 
dynamics codes are predicting a much higher edge velocity than is measured. This higher 
edge velocity on the suction side of the foil produces a lower pressure and therefore higher 
lift force. 
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IFLOW Measured MSES 



Figure 4-4: Predicted Verses Measured Stagnation Point 

3. Figure 4-4 shows the leading edge pressure distribution predicted by IFLOW and MSES 
as well as the experimentally measured leading edge data. It is clear from this figure that 
the observed stagnation point is more on the suction side of the foil than the predicted 
stagnation point. This conclusion is consistent with the fact that the experimentally 
measured lift is lower than predicted. 

4. Figure 4-5 shows the leading edge pressure distribution predicted by MSES when the angle 
of attack is adjusted to force the lift coefficient to agree with measurements. The updated 
stagnation point matches closely with the measured stagnation point as is expected. The 
angle of attack adjustment required to cause this change in lift is -1. degree. 




Figure 4-5: MSES Stagnation Point when Cl is Matched with Experimental Results 
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In summary, it is virtually certain that there are significant differences between predicted 
and measured flow velocities. The differences in lift are not a function of the method used to 



measure lift, but an actual difference in the flow field. 



4.4 Two Dimensional Flow 

Another potential difference between the flow codes and the experiment is a difference in aspect 
ratio. The computational fluid dynamic codes analyze a two-dimensional solution which uses 
an infinite aspect ratio assumption. The experiment uses actual foil aspect ratio of 0.9 and uses 
the tunnel walls to simulate an infinite aspect ratio. These sidewalls actually experience some 
boundary layer growth and therefore could lead to three dimensional effects in the tunnel. 

4.4.1 Gross Flow Field 

Jorde [l] measured the flow velocities across the tunnel and showed the flow at a given location 
in the gross flow field is two dimensional. 

4.4.2 Foil Boundary Layers 

Although the gross flow field is shown to be two dimensional, it is possible that there are three 
dimensional effects in the boundary layers which cause changes in boundary layer growth on 
the foil. These effects could result in significantly different lift, particularly considering that 
there is a separated flow region near the trailing edge of the foil. In order to eliminate the 
interaction between the tunnel wall boundary layers and the foil boundary layers fences are 
constructed out of aluminum and put on the foil to block any travel of the wall boundary 
layer down the span of the foil. This fence geometry is shown in Figure 4-6. When the actual 
wake velocities were measured with and without the fences there was no change. 

It is possible that the fences used here developed a boundary layer of their own, and then 
the fence boundary layers interacted with the foil boundary layers to cause a similar effect as 
the tunnel wall boundary layers. In order to test this hypothesis a set of modified fences were 
constructed that covered only the last 10% of the foil where there is significant thickness to the 
boundary layer. This also covers the region where separation is expected. Again the wake 
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Tunnel Wall 




Figure 4-6: Fence Geometry Used to Evaluate Sidewall Boundary Layer Effects 
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Figure 4-7: Modified Fence Geometry 
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Figure 4-8: Suction Side Tunnel Wall Boundary Layer Measured vs Predicted 

parameters and measured lift coefficient were unaffected. The conclusion is that there is no 
significant interaction between the wall boundary layers and the foil boundary layers. 

4.4.3 Tunnel Wall Boundary Layers 

It is conceivable that the tunnel wall top and bottom boundary layers would not behave as two 
dimensional walls due to interactions with the tunnel sidewall boundary layers. In particular, 
if there were significantly more growth on the suction side wall than predicted, it could result in 
an effective "‘decamber mg'’ of the foil. This would make the foil seem as if it had less camber 
than the actual geometry and reduce the experimental lift. 

This hypothesis is tested by measuring the tunnel suction side wall boundary layer at the 
front and back of the foil and comparing these measurements with the predictions of the same 
parameters in MSES. Figure 4-8 shows that there is good agreement between predicted wall 
boundary layer behavior and measured wall boundary layer behavior. 

4.4.4 Summary 

At this point in the analysis of the highly cupped trailing edge geometry the major experimental 
uncertainties have been addressed, but there is still a very large difference between measured 
and predicted performance. On the basis of this, it appears the numerical methods are not 
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satisfactory for this application. If the numerical methods are truly at fault, the errors here 
are clearly unacceptable. Given the presumption that two dimensional hydrodynamics can be 
accomplished by computational fluid dynamics regardless of geometry, there is still a hesitancy 
to trust these results. 
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Chapter 5 



Analysis of Blunt and Streamlined 
Foil Sections 



Based on the previous chapter, computational fluid dynamics and experimental results have 
important differences that result in drastically different lift predictions for the cupped foil 
geometry. The question that then arises is whether similar issues exist between CFD and 
experiments for other types of foil geometries. Experiments are conducted on a variety of foil 
shapes and computations are performed in MSES. 

5.1 Traditional Streamlined Sharp Trailing Edge 

The streamlined foil geometry features reduced loading near the leading edge for cavitation 
performance and a typical unloaded, streamlined and sharp trailing edge. This type of foil 
would be expected to be relatively easy for CFD codes to predict performance since the angle at 
which the flow leaves the foil is defined by the streamlined trailing edge. The experiments were 
conducted by Kimball [4] and the MSES computations were conducted for this thesis. The 
data taken by Kimball was analyzed using the algorithms applied on the rest of the foils in this 
experiment and the results were consistent with his conclusions. Figure 5-1 shows that there 
is almost perfect agreement between measured and predicted lift for this foil. The analysis 
conducted here validates the accuracy of MSES as well as the experimental setup for traditional 
foil designs. 
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Angle of Attack 



Figure 5-1: Numerical and Experimental Lift for Streamlined Sharp Trailing Edge Foil Geom- 
etry 



5.2 Blunt Trailing Edge 

5.2.1 Numerical and Experimental Results 

In the same experiment shown above Kimball conducted an analysis of a foil with a blunt 
trailing edge as shown in Figure 5-2. This foil was modeled after the parent foil discussed in 
the previous section, but includes a blunt trailing edge with a height of 2% of the chord length. 
The results for this foil are not as encouraging. There is significant error introduced due to 
the blunt trailing edge as shown in Figure 5-2. 

Another foil with only a 0.75% of the chord length blunt trailing edge is analyzed in MSES 
and experimentally. The differences between experimental and CFD predictions for this foil 
showed even greater error than the foil with more blunting at the trailing edge. The outcome 
of this analysis is shown in Figure 5-3. It is interesting to note that even though the extent of 
the blunting is smaller, the disagreement between experiment and CFD is larger. The nature 
of these results was somewhat surprising and caused the question to arise whether or not the 
boundary layer method of MSES would give the same results that would be obtained using a 
Reynolds Averaged Navier Stokes solver. 

DTNS2D is a generalized 2-D domain incompressible RANS fluid flow solver developed at 
the U.S. Naval Surface Warfare Center, Carderock Division. DTNS2D uses a finite volume 
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Figure 5-2: Numerical and Experimental Lift for Blunt Trailing Edge Foil Geometry 

finite difference formulation. It has been used previously for analyzing traditional surfaces with 
results that agree with experimental data [11]. The code DTNS2D was developed by Gorski 
[12]. The code utilizes a third-order upwind difference total- variation-diminishing (TVD) 
scheme applied to the convection terms and a second-order central difference scheme applied 
to the diffusion terms. The Baldwin-Lomax turbulence model and K-e models were both used 
for the analysis performed here and the results were consistent with 0.004 in lift coefficient. 

5.2.2 Analysis 

One major difference between the sharp streamlined trailing edge foil and the blunt trailing 
edge foil geometries is the separated region that occurs just beyond the trailing edge of a blunt 
foil. The obvious conclusion is that either the accuracy of the experiment, the accuracy of 
the flow code or the ability to represent the conditions in the flow code with an experiment is 
effected by this region. Here the latter of these possibilities is explored. 

It is possible that due to the change in geometry between the streamlined and blunt trailing 
edges that the ability to model the conditions used in CFD is lost. In order for this to be 
true there would need to be a significant global impact on the flow problem due to a very 
localized feature at the trailing edge. In the case of the cupped foil, this effect was eliminated 
by installing flow fences to block any communication between the foil boundary layer and the 
tunnel wall boundary layer. 
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Lift Coefficient 




Figure 5-3: Numerical and Experimental Lift for Modestly Blunt Trailing Edge Foil Geometry 
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modestly blunt foil C L = 0.457B 




Figure 5-4: Modestly Blunt Foil Boundary Layer Thickness 

In the case of the blunt foil, the sensitivity of foil performance to contamination of the 
boundary layer is studied and the contamination levels are measured. Figure 5-4 shows the 
predicted boundary layer thickness near the trailing edge from MSES and the experimentally 
measured boundary layer thicknesses. There seems to be good agreement between the experi- 
ment and predicted boundary layer thicknesses. 

In order to measure the sensitivity of foil performance to contamination of the boundary 
layer, false sources are placed in the separated region behind the foil. These sources are adjusted 
in magnitude until the lift coefficient realized in the experiment is matched. Another method of 
increasing the boundary layer thickness is reducing the Reynolds Number. When the Reynolds 
Number is reduced to 2*10 4 the lift coefficient predicted agrees with the experimental lift 
coefficient. Figure 5-5 shows the boundary layer displacement thickness predicted each way as 
well as the physically measured boundary layer thickness. It is clear that the measured value 
is closest to the MSES predicted boundary layer with no contamination at the appropriate 
Reynolds Number. This seems to indicate that although the predicted lift is wrong, the 
boundary layer calculations in MSES are functioning well. Based on this it would seem that 
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Figure 5-5: Wake Displacement Thicknesses Assuming Various Experimental Errors 
the tunnel is still providing a two dimensional environment for the foil to operate in. 

5.2.3 Effectiveness of Boundary Layer Trips 

Boundary layer trips are placed on the foil at 10% of the chord length to cause transition to a 
turbulent boundary layer. These trips are simulated in MSES by tripping the boundary layer at 
10% of the chord length. If these trips were ineffective it could conceivably effect the boundary 
layer size and therefore lift. This concept is explored on the modestly blunt trailing edge foil. 
Table 5.1 shows the variation in predicted lift coefficient for the conditions stated. It is also 
possible that the tunnel turbulence would force an almost immediate turbulent boundary layer 
which is modeled by putting the forced trip at the 1% position. These variations are small 
compared and generally in the wrong direction when compared to the experimental results. 
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Suction Side Trip 


Pressure Side Trip 


Lift Coefficient 


forced (0.1) 


forced (0.1) 


0.326 


natural (0.86) 


forced (0.1) 


0.409 


natural (0.86) 


natural (0.95) 


0.389 


forced (0.1) 


natural (0.95) 


0.301 


forced (0.01) 


forced (0.01) 


0.3244 


EXPERIMENT 


both forced (0.1) 


0.205 



Table 5.1: Lift Coefficient verses Turbulence Stimulator Effectiveness 
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Chapter 6 



Conclusions 



It has been shown in this thesis that experimental results for a two dimensional foil are accurate 
and repeatable. Many areas of uncertainty have been analyzed and none of the uncertainties 
can explain the inconsistencies between numerical and experimental foil performance. Based on 
these results, the following is a summary of the capabilities and limitations of CFD in predicting 
two dimensional foil performance: 



1. CFD is fully capabile of predicting performance of traditional foils with sharp unloaded 
trailing edges including cambered and symmetric designs. 



2. CFD is fully capabile of predicting performance of blunt trailing edge foils with no camber 
loading. 

3. CFD effectively predicts lift slope for asymmetric blunt trailing edge foils. 

4. ^FD over predicts the contribution of camber on asymmetric blunt trailing edge foils. 

5. 1 CFD over predicts the contribution of camber on foils with premature separation due 
to excessive trailing edge curvature. 



lr The offset in angle of attack is somewhat of a mystery, but a very significant effect. Continuing research will 
go into this area. The outstanding agreement between RANS codes and boundary layer codes seem to indicate 
that there is a common assumption between these programs that is in error. One current notion is that the 
organized asymmetric vortex shedding that goes along with a cambered separated trailing edge has a mean effect 
on the lift. If this were true, it may not be captured in a time average code. 
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6. 2 CFD over predicts the lift slope on foils with premature separation due to excessive 
trailing edge curvature. 



2 The inaccuracy in lift slope is most likely caused by incorrect prediction of the movement of the separation 
point as the angle of attack changes. 
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Appendix A 



Water Tunnel Three Dimensional 
Effects 



Figure (A-l) shows the effect of side wall boundary layer growth on control volume measure- 
ments. In a two dimensional control volume (shown in Figure (A-2)) the change in momentum 
between the fluid entering and leaving the control volume is measured at the boundaries to 
determine the net force on the control volume. These calculations utilize a unit thickness dz 
for determining the mass flow into and out of the control volume. If there is any net flow into 
or out of the control volume from the sides, this is not taken into account by this boundary 
calculation. 

A "typical” boundary layer growth of 5mm over the length of the control volume is consid- 
ered. The tunnel width is 500mm, so the boundary layer growth represents 2% of the total 
tunnel width. Based on the side wall boundary layer growth alone an increase in momentum 
of 2% is expected. The view of the control volume in Figure (A-2) is considered. This side of 
the control volume has a top and a bottom that follow a streamline in this plane. There is a 
net reduction in momentum of fluid in this plane based on the viscous friction with the surfaces 
of the foil. For a typical drag coefficient of Cq = 0.01 and a typical control volume height of 
^ of the height of the tunnel there is a momentum deficit caused by the foil of 1.5%. 

The velocity and quantity of flow into the sides of a control volume are difficult to predict 
and cause significant, sources of errors. It has been shown here that these sources of errors are 
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Side Tunnel Wall boundary laye^ 





Side Tunnel Wall 




Figure A-l: Depiction of Two Dimensional Control Volume on Edge 




Upper Tunnel Wall 






Lower Tunnel Wall 
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Figure A-2: Control Volume Considered for Analysis 
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on the same order of magnitude as the parameters being measured. This explains why it has 
been the consistent experience at MIT that the control volume measurement of drag forces is 
unreliable and sometimes even shows a net thrust! 
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Appendix B 



Laser Doppler Velocimetry Data 




Lower Tunnel Wall 



Figure B-l: Data Summary Key 



Figure B-l shows the geometry of the control volume data plotted. Each of the plots 
(Figures (B-2), (B-3) and (B-4)) shows the properties indicated verses path length along the 
control volume. The data is plotted counter clockwise starting at the lower right hand corner 
of each control volume. 
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Normalized Static Pressure Normalized Velocity 



1.25 




' Normalized Horizontal Velocity 
Normalized Vertical Velocity +1 




Normalized Static Pressure 

Figure B-2: Control Volume #1 Raw Data 
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Normalized Static Pressure Normalized Velocity 
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Normalized Static Pressure 

Figure B-3: Control Volume #2 Raw Data 
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Normalized Static Pressure Normalized Velocity 



1.25 



1.13 



1 



0.88 



0.75 

0.63 



0.5 

0 0.8 1.6 2.4 3.2 4 

Path Position (S/C) 

Normalized Horizontal Velocity 
Normalized Vertical Velocity +1 





Normalized Static Pressure 



Figure B-4: Control Volume #3 Raw Data 



56 



Appendix C 



PIV Experimental Data 



57 




Figure C-l: Particle Image Velocimetry Experimental Setup 

Particle image velocimetry is an experimental method which consists of lighting a plane of 
fluid containing flourescent particles with a laser sheet and taking pictures of their position 
with a camera. Two pictures are taken close in time with one another and the distance of the 
movement of the particles between laser pulses combined with a known time between laser pulses 
gives a snapshot of the velocity at an instant in time. For this experiment the time between laser 
pulses was 0.037 milliseconds and the camera area covered was about 5 centimeters. Figure 
C-l shows the experimental setup used for these PIV measurements. Measurements are only 
taken on the suction side of the foil in this portion of the experiment. 

The particle image velocimetry measurements at each instant in time are averaged to provide 
a mean velocity. In this experiment only 126 individual velocity values were used to determine 
the average. It is clear from laser doppler velocimetry experience that 126 samples in an average 
is not sufficient to provide a reliable mean value in turbulent areas. Another weakness of the 
particle image velocimetry measurements in this experiment is that each sample was taken in a 
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Figure C-2: Summary of PIV Data Including Pressure Contours 

single burst of laser pulses and pictures. This means that any tunnel speed fluctuations at low 
frequency (on the order of seconds) will show up in the particle image velocimetry measurements 
as a shift in the mean values for the particular sample. In short, the mean values of the particle 
image velocimetry measurements in this experiment are suspect, but the trends shown in each 
sample are likely to be reliable. 

Figure C-2 shows the datafield considered in the particle image velocimetry portion of 
the experiment. The contours shown are pressure contours. The boundaries shown are the 
individual pictures that were taken of the different foil sections. Discontinuities from one picture 
to the next are caused by the variations mentioned above. Each of these pictures contains a 
128 by 128 grid of horizontal and vertical velocity at each particular location considered. 

The following diagrams are plots of the boundary layer parameters from the various sources 
that are available. The PIV boundary layers have been scaled to match the projected LDV 
measurements since the mean value of the LDV measurements are highly reliable. These 
diagrams are discussed further in the body of this thesis. 
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